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Abstract 


Fuzzy multiobjective linear bilevel programming (FMOLBP) problems are 
studied in this paper. The existing methods replace one or some determin- 
istic model(s) instead of the problem and solve the model(s). Doing this 
work, we lose much information about the compromise decision, and it does 
not make sense for the uncertain conditions. To overcome the difficulties, 
Zadeh’s extension principle is applied to solve the FMOLBP problems. Two 
crisp multiobjective linear three-level programming problems are proposed 
to find the lower and upper bound of its objective values in different levels. 
The problems are reduced to some linear optimization problems using one 
of the scalarization approaches, called the weighting method, the dual the- 
ory, and the vertex enumeration method. The lower and upper bounds are 
estimated by the resolution of the corresponding linear optimization prob- 
lems. Hence, the membership functions of compromise objective values are 
produced, which is the main contribution of this paper. This technique is 
applied for the problem for the first time. This method applies all informa- 
tion of a fuzzy number and does not estimate it by a crisp number. Hence, 
the compromise decision resulted from the proposed method is consistent 
with reality. This point can minimize the gap between theory and practice. 
The results are compared with the results of existing approaches. It shows 
the efficiency of the proposed approach. 
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1 Introduction 


A bilevel multiobjective programming problem (BL-MOPP) is applied to 
model decentralized decision problems that contain the objectives at the 
upper level, called leader, and the objectives at the lower level, called the 
follower. The decisions are made in a hierarchical order from the upper 
level to the lower level. The decision-makers in each level try to optimize 
their objective functions. However, the decision in a level may be affected 
by the objective values in another level. The objectives in both levels may 
have confliction objectives, which should be optimized simultaneously. The 
coefficients in the objective functions and constraints are determined by ex- 
perts; see [3, 46]. In the most real-world situations, the possible values of the 
coefficients are usually estimated imprecisely or vaguely by experts. These 
values cannot be presented in terms of crisp values because we miss much 
information. The most appropriate form to express the values is fuzzy sets. 

Many researchers have been studied bilevel programming problems since 
it was introduced by Von Stackelberg [36]. A class of algorithms and ap- 
proaches was well designed to solve the problem like the Kuhn—Tucker method 
(7, 8, 9, 16], the vertex enumeration method [9], the K-best approach [9, 2], 
and the branch and bound procedure [14]. The BL-MOPPs were studied to 
formulate the problems of the real-world, for example, the location planning 
problem for stone industrial parks; see [5]. Their coefficients are expressed 
by crisp numbers. On the other hand, the coefficients cannot usually be 
expressed in terms of a real number in the uncertain conditions. Also, mul- 
tilevel programming problems (MLPPs) in [15, 35, 34] were discussed for 
hierarchical decentralized planning problems. An overview of bilevel pro- 
gramming problem was provided by Colson, Marcotte, and Savard [13]. The 
fuzzy approach for MLPP in [34] was extended to solve bilevel and three-level 
nonlinear multiobjective programming problems by Abo-Sinna and Baky [2]. 
Interactive fuzzy programming was developed to solve fuzzy multilevel linear 
programming problem in [29]. The balance space method was proposed for 
the nonlinear multiobjective bilevel programming problem in [2]. The fuzzy 
goal programming (FGP) algorithm was used to solve BL-MOPP in [6]. An 
algorithm was designed based on FGP to solve MLPP in [9]. An FGP model 
was proposed to solve BL-MOPP by Pramanik and Dey [4]. This model does 
not follow the hierarchical structure of bilevel programming and does not 
consider the upper level decisions, and it is as a single-level multiobjective 
programming problem. In bilevel programming problems, the decisions of 
leader should be dominant over that of followers. Mohamed [5] well devel- 
oped an efficient method to solve the multiobjective programming problems 
using FGP approach. Moitra and Pal [15] then extended the concept of FGP 
approach to finding a satisfactory solution for bilevel programming problems. 
Moreover, FGP approach was extended to solve a deterministic decentralized 
BL-MOPP in [6] and MLPP in [4]. Due to the existence of uncertainty in 
real-world situations, some or all parameters in objective functions and/or in 
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constraints of the leader and/or the follower can be described by fuzzy num- 
bers. Budnitzkia [5] designed a solution procedure to solve the linear fuzzy 
bilevel programming using a-cuts of the fuzzy polytopes. Zhang and his 
colleagues [17, 35, 37] utilized the extended solution concept and theorems 
of bilevel programming and presented procedures to solve the fuzzy linear 
bilevel programming (FLBP) in special cases of membership functions such 
as the triangular and the general form [42, 43, 44, 45]. When the leader, the 
follower, or both of them have multiple objectives, the fuzzy set technique 
was applied to investigate FLBP problems by Zhang, Lu, and Dillon [46]. 
They also extended their previous research to the fuzzy multiobjective linear 
bilevel programming (FMOLBP) problem and developed an approximation 
branch-and-bound algorithm to solve the FMOLBP problem; sce [46]. A new 
definition proposed in [32, 35] implies that the upper level constraints con- 
taining the second variable move to the lower level. This point changes the 
nature of the problem; see [1, 8]. On the other hand, the method given in 
[46] solves several crisp multiobjective nonlinear programming problems along 
with many crisp objectives and constraints by the approximation branch and 
bound method until the termination criterion holds. Hence, it will have high 
computational complexity and we have to check the complementary slackness 
conditions in each of the iterations, which is time-consuming. 


Toksari and Bilim [13] focused on decentralized bilevel multiobjective frac- 
tional programming problems (DBL-MOFPPs) with a single decision-maker 
at the first level and multiple decision-makers at the second level, where the 
coefficients are crisp numbers. They presented an FGP based on Jacobian 
matrix for DBL-MOFPP. In the approach, membership functions of fuzzy 
goals are constructed for all objectives at two levels, and they are linearized 
by using the Jacobian matrix. Then the FGP approach is applied to ob- 
tain the highest degree of each of the membership goals by finding the most 
satisfactory solution for all decision-makers. This algorithm was applied for 
the crisp problems based on FGP. Peric, Babic, and Omerovi [15] studied a 
similar problem with crisp coefficients and presented a method based on the 
FGP. The approaches fail to solve the proposed model of this paper due to its 
fuzzy coefficients. Liu and Yang [18] proposed an interactive programming 
method to obtain the compromise optimal solution of the multilevel multi- 
objective linear programming problem with crisp parameters. The resolution 
process is done in two stages: analysis and decision-making. If the objective 
values were not satisfactory, then the decision-maker improves the values by 
giving a concession to the unsatisfied values. When the satisfactory degree 
in upper levels is mat, the problem of lower levels will be solved. The pro- 
posed method in [18] was designed for the case that the coefficients are crisp 
numbers. The proposed method cannot be applied to the proposed model 
with fuzzy coefficients. Kamal et al. [19] applied crisp bilevel multiobjective 
programming for production planning problems. Then an FGP algorithm 
was designed to solve the problem. This approach cannot also be applied 
to solve FMOLBP due to fuzziness of its coefficients. Abdelaziz and Mejri 
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[23] formulated a shared inventory model as a bilevel programming problem. 
Both emergency and backorders are considered in the model. The interaction 
between decision-makers is done in an uncertain environment of probability 
type. They designed a decentralized bilevel programming problem, where 
the leader has multiple objectives. The numerical study was focused on the 
simulation. We cannot use the approach to solve the proposed model with 
fuzzy coefficients because the approach in [23] discusses with the uncertainty 
of the probability type. Baky Eid, and El Sayed [3] extended the FGP to 
solve BL-MOPP with fuzzy demands that are given as triangular fuzzy num- 
bers. The FGP algorithm is applied to attain the highest degree of each of 
the membership goals by minimizing their deviational variables and does not 
present the membership functions of the objectives of leader and follower, 
completely. In this paper, we will answer the following questions about the 
proposed method in [3]: 


1- How can we find the membership functions of objectives of leader and 
follower instead of obtaining the highest degree of each of the membership 
goals? 


2- Finding the required goals is not always easy, especially, when the di- 
mension of the problem increases. How can we overcome these difficulties? 


The answers to the above questions are the motivations of this paper. In 
this paper, we will consider the FMOLBP problem in a general case, where 
all coefficients are fuzzy numbers. Zadeh’s extension principle [6, 40, 41] is 
applied to solve the problem. In this approach, a pair of crisp multiobjective 
linear three-level mathematical programming problem is formulated to com- 
pute the lower and upper bounds of the a-level of the objective values. Then, 
one of the scalarization methods, called weighting method, dual theory, and 
vertex enumeration approach are applied to solve the crisp multiobjective 
linear three-level programming problems. Finally, the membership function 
of the fuzzy compromise objective value is derived numerically by enumer- 
ating different values of a. This is the first contribution of the paper. In 
this approach, the KKT conditions are not used because of its nonlinearity 
nature. When the constraint functions at the upper level have an arbitrary 
linear form, the K-best algorithm cannot always find the optimal solution; 
see [31]. Hence, the most suitable method will be the vertex enumeration ap- 
proach along with the weighting method and dual theory. These techniques 
preserve the linearity property of the problem. In addition, the results of the 
proposed approach are compared with the results of existing approaches. It 
shows the efficiency of the proposed approach with respect to the existing 
methods. This point is the second contribution of the paper. 

The organization of this paper is as follows: Section 2 is formed of three 
subsections. The notations and properties of fuzzy numbers are reminded. 
Also, the weighting method is reviewed to solve multiobjective linear pro- 
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gramming problem in the second subsection. In the third subsection, the 
multiobjective linear bilevel programming problem is formulated and an ap- 
proach is presented to solve it based on the weighting method, dual theory, 
and vertex enumeration approach. Section 3 is divided into two subsections. 
The first subsection formulates the FMOLBP problem. The second sub- 
section designs a solution procedure to solve the problem based on Zadeh’s 
extension principle. Section 4 designs an efficient algorithm to solve the prob- 
lem of FMOLBP based on the results of Section 3. A numerical example is 
presented to illustrate the algorithm in Section 5. Section 6 presents a com- 
parison among the proposed algorithm and the existing methods to show the 
efficiency of the proposed algorithm. Conclusions are presented in Section 7. 


2 Preliminaries 


This section provides the required preliminaries in terms of three subsections. 
The preliminaries are about fuzzy numbers, weighting method in multiobjec- 
tive linear programming problem, and multiobjective linear bilevel program- 
ming. 


2.1 Fuzzy numbers 


Let R be the set of all real numbers, let R” be n-dimensional Euclidean 


space, and let x = (21,...,%n),y = (y1,---, Yn) € R” be two vectors, where 
vi,yi € R, i =1,...,n. For any two vectors x,y € R”, we write x > y if 
and only if x; > y, for alli =1,...,n. The notation “o” displays the zero 


vector. Its dimension is determined in the context. 


Assume X C R. In this paper, it is supposed that the set X is equal to R 
or [0,1]. If @ is a fuzzy subset of X with membership function pg : X > [0,1], 
then the a-level set of @ is as Gg = {a € X | fa(x) > a} for all a € (0, 1] and 
ag = cl({x € X | a(x) > Of). 


Definition 1. (a) Assume that FN(X) is the set of all fuzzy subsets @ of X 
with pg satisfying the following conditions: 

(i) There exists x € X such that a(x) = 1; 

(ii) wa(A.x + (1 — A).y) > min{pa(x), ua(y)} for all vy © X and 2 € (0, 1]; 
(iii) {@ € X | wa(x) > a} = Ga is a closed subset for each a € (0, 1]; 

(iv) The 0-level set Go is a compact subset of X. 

(b) A triangular fuzzy number is a fuzzy number that is presented by @ = 
(a1, 42,43), where a, < ag < ag and its membership function is as follows: 
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r—ai 4 
=. ifa, <a@ < ag, 
= — a3—& } 
ba(x) = =o if ag <a < as, 
0 otherwise. 


With regard to Definition 1, the a-level set a@ can be written as dg = 
[aX a]; see [40]. In the triangular case, the a-level set @ is as Ga = [a1 + 
(a2 — a1).a, a3 — (a3 — a2).a]. If @ is a crisp number with the value m, then 


its membership function is as follows: 


mats f ife=m, (1) 


0 otherwise. 
The set of all finite fuzzy numbers on R is denoted by FN*(R). 


Definition 2. Let a; ¢ FN(R) for i=1,...,n. We define @ = (@,...,@n) 
and pig : R" > [0,1] by wa(x) = min pa, (xi), where x = (21,...,2%,) € R”, 


peeey 


and @ is called an n-dimensional fuzzy number on R”. If a; € FN*(R), 
a= 1,...,n, then a is called an n-dimensional finite fuzzy number on R”. 
Let F.N(R") and F.N*(R”) be the set of all n-dimensional fuzzy numbers 
and the set of all n-dimensional finite fuzzy numbers on R”, respectively. 


2.2 Weighting method in multiobjective linear 
programming (MOLP) problem 


Consider a multiobjective linear programming (MOLP) problem 


min 2(#) = (z1(2),..., 2n(2)), (2) 
st. cEX CR", (3) 


where x € R” and X denote the decision-making variables and feasible set, 
respectively. Moreover, Z : R" —+ R* and z;(x),i=1,...,k, are real-valued 
functions over X. 


Definition 3 ({17]). The vector x* is said to be a Pareto optimal solution 
(or compromise solution) if and only if there does not exist another x € X 
such that zi(x) < z(x*) for alli and z;(x) 4 z;(a*) for at least one 7. 


The weighting method for obtaining a Pareto optimal solution is to solve 
the weighting problem formulated by taking the weighted sum of all of the ob- 
jective functions of the original MOLP (1)-(3). Thus, the weighting problem 
is defined by 
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min W.z(#) = Yom zi(x (4) 
k 
st. «vex, yet w,>0, for i=1,...,k, (5) 
i=l 
where W = (wi,..., Wx) is the vector of weighting coefficients assigned to 
the objective functions, and it is assumed that W = (w1,...,wz) > 0. The 


following theorem presents a relationship between the Pareto optimal solution 
of problem (1)—(3) and the optimal solution of problem (2)—(5). 


Theorem 1 ({17]). If* € X is an optimal solution of the weighting problem 
(2)-(5) for some W > o, then x* is a Pareto optimal solution of the MOLP 
(1)-(3). Conversely, if 2* € X is a Pareto optimal solution of the MOLP 
(1)-(3), then z* is an optimal solution of the weighting problem (2)—(5) for 
some W = (wi,...,Wk) > 0. 


2.3 Multiobjective linear bilevel programming(MOLBP) 


First of all, throughout this paper, we assume that X = {a € R"|x > o} and 
Y = {y € R™|y > 0}. The MOLBP problem is formulated as follows: 
Forzé€ X CR", yEY CR", F:X xY —R‘,and f:X x Y — R', we 
define 


/ / / Uy 
Z(c1,.-- »Cs,d1,- an ,ds,C4,- +) Cz, y,.. .,d;,A1, By, Ao, Ba, k1, ke) 


= min F(z, y) =(c,.4+d1.-y,...,€5.2+ds-y)?, (6) 
st. Ay.at+ Bry < ky, (7) 
min f(x, ¥) =(c.otdy,...,c.2+d,y)*, (8) 
s.t. Ag.a+ Boy < ke, (9) 
where 
C= Gases Gu) G = (Gye) OR tH] = lest 
= (dyi,...,dma), & = (djs---5 Uj) ER”, §=1,---,8, fale t 
ki = (ki1,..-,kp1) € R?, ko = (Kia,...,kg2) € RY, Ar = [aij] € R°*”, 


B= [bi, i] € € Re*™, Ag = [ei] E Res", and Bj = [si9] € RIX™, 

Let S = {(a,y) € X x Y|A1.a + Bry < ky, Ao.x + Boy < ke} and 
S(a) = {y € Y | Boy < kp — Ao.a}. Denote by R(x) the set of the Pareto 
optimal solutions of the lower level problem for any fixed x. Then we have 
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the following definitions. 


Definition 4 ({8]). For a fixed (x,y) € S, if y is a Pareto optimal solution to 
the lower level problem, then (x, y) is a feasible solution to problem (6)—(6). 


Definition 5 ((8]). If (x,y) is a feasible solution to problem (6)—(6) and z is 
a Pareto optimal solution to the upper level problem for fixed y, then (a, y) 
is a Pareto optimal solution to problem (6)—(6). 


We now apply the weighting method and rewrite problem (6)—(6) as fol- 
lows: 


TF / / Uy Ul 
Z(C1,---5€s,d1,---,ds,C,,---,G,d4,---,d;, A1, Bi, Ao, Bo, ki, ke) 


=min WF = i (cj. + diy), 10 
min (x,y) 2 (ci. + diy) (10) 
s.t. Ay2+ Bry Ski, Sw =1,w; > 0, for t= 1.06458) (11) 
i=1 
t 
in W’. — (C.2 + dy. 12 
min W".f(x,y) Di (Gj.a + d;.y), (12) 


t 
s.t. Ao. + Boy < ko, YS) wi =1,w; > 0, for i=1,...,t,. (13) 
i=1 


For some fixed x > 0, the compromise solution of the lower programming can 
be found by solving the following linear programming with a given weight 
W'>o. 


t 
min we wy.(c,.2 + di.y), (14) 
y ; 
t=1 
s.t. Bo.y < ko —= Ao.2, (15) 
t 
Swi = 1, uj > 0, for i=1,...,t,y>o. (16) 
i=1 
Since cj.2, fori = 1,...,¢, are constant, we can ignore the terms. Hence, a 


dual multiplier wu; is considered for each constraint (15), the primal variables 
are in terms of the vector of y, and other expressions are considered as con- 
stant. Thus, we can get the dual of problem (14)—(16) according to the dual 
definition for linear programming problems in [40] as follows: 
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max (kz — Ag.a)? .u, (17) 
t 
s.t. BE wu > So wi.d, (18) 
i=1 
fl 
So =, w,>0, for i=1,...,t, u>o, (19) 
i=1 


where u € R¥® is the dual multipliers. The following theorem converts 
MOLBP problem (6)-(6) to a one-level programming problem with one ob- 
jective function. 


Theorem 2. The pair (2*,y*) is the compromise solution of problem (6)— 
(6) if and only if there exists u* such that (2*,y*,u*) is the solution of the 
following program: 


min F(z,y)= S- w;.(c.c + di-y), (20) 
i=1 
S.t. Aj .x + By.y < ky, (21) 
t 
By.u> > widi, (22) 
i=1 
t 
S¢ wi.(dj.y) — (ke — Ao.2)".u =0, (23) 
i=1 
t 
Sui = 1, uj > 0, for i=1,...,t, Zy,u>o. (24) 
i=1 


Proof. The proof is similar to the proof of Theorem 1 in [33]. 


Let U = {u|B2.u > sae w}.di.,u > o} denote the feasible region to linear 
programming (17)—(19). If U #0, then the set U has at least one vertex and 
at most finite vertices. Moreover, if there is an optimal solution to the linear 
programming (17)—(19), then it must be one of the vertex points of set U (see 
[28]). With attention to the above results, we can convert problem (20)—(24) 
to a series of the following linear programming problems by obtaining all 
vertices of U, denoted by U” = {u',u?,...,u"}, according to the method in 
linear programming problem from [28]. 
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LP(u?): min F(a,y) = Yel c.« + dj-y), (25) 
xy 

S.t. Ay.“ + ae < ky, (26) 

Yow (d',.y) — (ky — Ag.x)? u? = 0, (27) 


Seale for 1=1,...,t, wy >o. (28) 
i=1 


We can solve the above LP(u?) forp =1,...,r. Let J C {1,...,r} denote 
the Indices of problems that LP(u'), for i € I, has an optimal solution. For 
i € I, let (x*, y") be the optimal solution for problem LP(u') and F(x", y*) = 
min{F(a',y’)|i € I}. Therefore, we have the following theorem. 


Theorem 3. The pair (x*,y*) is an optimal solution for problem (20)—(24) 
and a Pareto optimal solution for problem (6)—(6). 


We are now ready to formulate and solve fuzzy multiobjective linear 
bilevel programming. 
3 Fuzzy multiobjective linear bilevel programming 


In this section, the FMOLBP problem is formulated and a procedure is de- 
signed to solve the problem based on Zadeh’s extension principle [6, 40, 41]. 


3.1 The formulation of FMOLBP 


For € X CR",yeY CR”, F: X x Y — FN*(R®), and f:X x Y 
FN*(R‘), the problem of FMOLBP is formulated by 


Z= min F(a,y) = (@.¢+dzy,...,é.0 +ds.y)’, (29) 
rE 
s.t. Ay.2 + By.y < ka, (30) 


min f(r,y) =(Gat+dyy,...,g.0+d.y)’, (31) 
y 


S.t. Ag.x + Bo.y < ko, (32) 
E20, y>o, (33) 


= ins Oal es = lignite Ph ta haan 7 Sha gh 
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C= ise = yee dg EE RH is Set 
ky = (ki1,.--, kp) € FN*(R?), ko = (kiz,..., ha) € FN*(R), 
A = [aijlpxn, dig € FN*(R), Bi = [0; eee: ij € FN*(R), 


Ao = [Eijlaxn, €i3 € FN*(R), a By = [Sijlaxm, 513 € FN*(R). 
The feasible solution set of FMOLBP is defined as follows: 


S= {(z,y) € X x Y|Ay.2 + Bry < ky, Ag. + Bo.y < ko}. (34) 


In an FMOLBP problem, for each (x,y) € S, the value of the objective func- 
tions F(2,y) = (Fi(@,y),-..,Fe(x,y)) and fle,y) = (file,y), ++ fel@.y)) 
of leader and follower is s-dimension and t-dimension fuzzy vectors, respec- 
tively. 


3.2 The solution procedure for the problem (29)—(33) 


We will now consider problem (29)—(33). Let 

Marsrs +s Hepes Meijy ++ Megs Mayr Madge Mates Ha b= 1. - 58 7 = 
ee 
My t= 1,.--57RI = 1,-. ,m, fe,;, 7 = 1,...,¢,9 = 1,...,n, and pz,,, 


7 = 1,. et j = 1,...,m, denote the membership factions Cliy+ ++ Cniy 
Bisel a DM eswhanitlege nr ae oe eee sae eee ae eee 
on ee i ee ee 
1,...,m, ij, i= 1,...,4,9 =1,...,n, and 53,7 =1,...,¢,7 = 1,...,m, 
respectively. 


We have 0 = {(v,o(v)) | v € S(v)}, where 0 is one of the parameters 
Clis+++5Eniy Crayons G ng? na 7 nea: ae 4 = Les 505 By => Liveek by 
ka, kya, = Li iceg Dg = => Woe seg Gis, i= We cieg D594 => Me due IV bij, t= 
1,...,2,9 =1,...,m, @ij, ¢ = 1,...,¢,79 = 1,...,n, 67,1 = 1,...,4,9 = 
1,...,m, and $(@) is the support of 6. Denote the a-cuts of 0 by 


()a = [(8)q, (8)a] = [min{u € S(@) | we(v) = a}, 
max{v € S() | wo(v) > a}], — (35) 


where v is one of the above fuzzy parameters. These intervals show that the 
coefficients of objective function and constraints in problem (29)—(33) lie at 
the a-level set. Here, we intend to obtain the membership function Z. To do 
this, Zadeh’s extension principle [6, 40, 41] is used. Based on the extension 
principle, wz can be written by 
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z(z)= sup — mint He, (Cri), ++ Ming (Cri), Hay, (Cig)s +++ Ber, (Cry) 
ej,0%,d;.d', 
ky, Ro, AL 
vine ae 


May, (Ari)s +++ 1a, (ini), May, dig)s +++ May, (ding)s Mie, (Rit); 
Hijo (kj2), ba; (aij), HG, (dij), Me, (e4j), Ma,; (8:3), for all 1,J | 
Z£>= Z(c1,-- ‘Cay ets 24 day Chg+x' Ch, nad ..,d, Ai, By, Az, Bo, ki, ke)}, 


(36) 
where matrices c; = Se C5 = le ilikns d= du liens Oy = [d plixca; 
i= hiss 337 = = Li oa ky = = [kit |px1; kg = = [kj2|qx1; A, = = [aijlpxn; 
By = bis levees Ag = elon and Bz = [sij]qxm, are the supports of 
Cin Oy di, d';, $= en 83-7 = Thee: t, ki, ke, Ay, Ba, Ao; and Bo, respectively. 
Also ne: Nin: a ame y+) dh, Ar, Bi, Ao, Bo, ki, kz) is de- 


fined as follows: 


BCigcen ag Oey Gigs on day Cy pos 09s toe++>d,,A1, Bi, Ao, Ba, ki, kz) 
= min F(a, y) = (c..4 + dy.y,...,¢s.c+ds-y)’, 

rTExX 
s.t. Ay.2+ By.y < ky, 


min f(x,y) = (c.c+dy.y,...,¢.0 +d,y)*, 
yeEY 


If the a-cuts of Z at all a values are converted to the same point, then the ob- 
jective functions are crisp numbers. Otherwise, they are fuzzy numbers. Ac- 
cording to (36), 4 is the minimum of p45 (v) for v = Ci: Ci di, di, kin, kyo, aiz, biz, 
€ij, 8ij, for all 7,7. We need Ls (v) > a, for v = ci, Ch ds, dj, kit, kj2, Giz, bij, eij, 
sij, for all i,j, and at least one s(v) > a, for v = ci, cj, di, dj, kin, ky2, aig, bij, 
7 sij, for all 2,7, equal to a such that z = Z(c,... <oueilis recta Oa oadinsers 

,d,,A1,.Bi, Ao, Bo, ki, ke) to satisfy wz(z) = a. To compute the mem- 
ene function 7, we intend to compute the left and right shape function 
of wz. It is equivalent to finding we and ZY of the ie cuts of Z. Since ZL is i 
the minimum of 21; 2025 Cg Qj. 05 Mas Os onc9 ep days di, A, By, Ap, Bo, 
ki, ky) and ZY is the maximum of Z(c1,... sacle a cae aceig ey Qayy crepes 
Ai, By, Ao, Bo, ky, ka), we have 


Sf ri ! 
Ze =min{Z(c1,...,¢s,d1,...,ds,Cy,--- es oe .,d,, At, By, Aa, Bo, k1, ke) 
~\L ~\U ! 
| (a) <u< i i ire Neem 
7 Ul «+ 
1jr1+ + Amys Gigs dij, Cig, Siz, Kir, Kyo, for all i,j} (38) 


and 
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7U i too / 
Lo = max{Z(c1,...,¢s,d1,---,ds,Cy,+--5Gs ty-+-34;, At, By, Ag, Bo, ki, ke) 


“af; ~U U 
[ta SOL, Where. WU = Cie oes Ge igsess pleas Giger oss lays 
Ljr+++s dmg Aig, big, Cay, 813, Kir, kyo, for all 4, 7}. (39) 
They can equivalently be rewritten as two problems below: 
y y 
: d eae 
min(c,x =F 1Y, +++, €g& = sY) ’ 
cEeX 
SL a s.t. Aix “fe By < ky, 
Zy = min : ! / ! 1 \T 
(EL <v < (YU, where min(cja + djy,...,q¢@+dyy) , 
U = Class sCnir dias: + dma yeY 
Cha yeeey el ogaitlhas en a idly 
ere, 2a jay for all i,j. st. Aga + Boy < ka, 
(40) 
and 
d dsy)" 
min(cjx + diy,...,¢Cst +dsy)’ , 
rex 
SU s.t. Aix + Byy < ky, 
Za = max : ! 1 ! i \T 
(@E <u < (UY, where min(c}2z +diy,...,q¢@2+d,y)’, 
US CZs ion 4 Cpgs Olga ee 4 yds yEeY 
e415 Sa eng Ay daonak qing? 
aij. bij, e471 8ij>ki1,k7j2, for all i, j. s.t. Aga + Boy < ko. 
(41) 


We can use the weighting method for the problem to find the compromise 
solutions of problem (40). Hence, we can rewrite problem (40) as follows: 


S.t. Ax + Byy < ky, 
s 


Se is for t= Ly cg 


i=1 
L. = max t 
(@)E <u < (#)U, where . pes Vea 
= Lin Omir Airs 5 dmis man w;(c,-@) + 5 w;(d;-¥); 
Chgreees eng Agrees qaj? yeY fay ‘= 
ajgsbijseigr Sigs hits kjQ, for all i,j. = = 


S.t. Aox + Boy < ko, 


t 
Sa = he 0 fort = Tg 
i=1 


(42) 
Since model (42) finds the minimum of all minimum objective values, we 
should consider the objective function of leader in the best case, that is, we 
set c; and d; to their lower bound, in other words, (é;)% and (d;)z, for each 
i, respectively. Therefore, we can rewrite problem (42) as follows: 
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s s 
mip) wi((&)g-2) + d. wi((di)e-y); 
s.t. Ajx+ Byy < ki, 
8 


eS 1G 0, fort 1.48, 


eo i=1 

Zw = max t t 
(OE <u < (YU, where P ppt ish 
er ae min w;(c,.@) + ) w,(d;.y), 
045 bi3, e457 8ij> bis bjQ, for all 4, J. ye ae ar 


s.t. Agr + Boy < ka, 

t 

Yo SS 0, fort 1a ts 
i=l 


(43) 


Notation 1. Assume that A = [4ij;]mxn is a fuzzy matrix with fuzzy com- 
ponents. Then (A)% = [(@i;)E]mxn and (A)Y = [(4ij)Z]mxn- 


a 


To obtain the minimum objective value, we should consider the values 
of Ai, By, Ao, Bo, ky, and ko in (Aq)? (Bie, (Ao)¥, (By)¥, (kz, and (k2)z, 
respectively, because we can obtain the largest feasible solution set, in this 
case. Hence, the least value is found for the objective function. To this end, 
problem (43) is converted to the following problem. 


s 


s 
: ~\L 5 \L 
mon wi((Gi)q-%) + Dd, wil(diard) 

i= i= 
at. (AiG e+ (Bi)oy Ss (ee, 

s 

Sou; =1,w; > 0, fori=1,...,s, 

Zy= min Sy 
@E<v< @)E where 


t 

j ' min S wi(c.a) + ‘(d,.y) 

V=Ch pres eng gree Qing? Pa oanG W\G;-Y), 
i=1 


for all j. yey rr 


s.t. (Ag)Ya + (Bo) y < (ke), 


t 
og S10, > 0 fort aaa, 
i=1 


(44) 


With regard to Theorem 2, we can write the inner problem as follows: 
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i=l i=l 
st. (Ai)2a+ (Bi)oy < (kg, 


s 


YG SH 1 SU, for tS 1) 1.48, 


t=1 
5h Fi 
a i 5 ~ 
weer sok, do wildi.y) — ((ka)e — (Anan)? .u = 0, 
where v=ch.,..., cy — 
dh grees dhs: for an ve t 
((Bo)8)7.u > So wi.dt, 
; i=1 
Salata thfatata.st 
w=1 
(45) 
Or equivalently, we can rewrite problem (45) as follows: 
<L 2 ud . 
Z, =min ))wi((G)d-2) +S wi((die-y), (46) 
i=1 i=1 
st. (Ay)Pat (Bi)by < (hg, (47) 
t 
do wi (di-y) — ((he)& — (Aa)ax)” uw = 0, (48) 
i=1 
- t 
((B2)0)".u > So) wi.di, (49) 
i=1 
(dio <d,<(d)¥, for alli=1,...,t, (50) 


8 t 
So wi =1, 5) uj = 1,04, 0) > 0, for i=]1,...,8, j=1,...,¢, 
i=l i=l 


L,Y, U> o. (51) 


We can now remove the constraint (50) and equivalently rewrite problem 
(46)—(51) as follows: 
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Z. = min) wi((@)e.0) + So wil(dén), (52) 
st. (Ai)Uat (Bly < (hz, (53) 


wiley) S ((ka)e — (Aa)oa)?.us Sl wi((@iey), (54) 


s t 
So wi =1,) > uj = 1,4, 05 > 0, for pa Ne 38s FS deg b, 
i=l i=1 


L,Y,uU> o. (56) 
t 
Let U = {u>o0| ((B2)¥)? u> S w.(d,)U}. Assume that the vertex points 
i=1 
of U are denoted by the set as UY = {u!,...,u"}. Hence, we should solve at 
most r linear programming problems as follows: 


LP(u*): 25 = min 7 wi((&)Ec) + Yo wil(die-n), (57) 
st. (Ai)Sa+ (Bi)ly < (hi)k, (58) 


wildy) < ((ke)% — (A2)f a)? .u* < wild ay); 


t=1 
(59) 
s t 
So wi = 1,50 w; = 1, w;, w; > 0, 
i=1 i=1 
for i=1,...,8, =1,...,t, zy>o, (60) 
where k = 1,...,r. Solving r linear programming problems LP(u"), we can 


find the optimal solution of the problem for each a € [0,1] as it was mentioned 
in Subsection 2.3. Finding the optimal values of variables x and y for each 
a € [0,1], we can estimate values (é;)4.2 + (d;)4.y for each a € [0,1] and 
i=1,...,s. Hence, we can obtain A 

We now focus on obtaining the ae. To do this, the weighting method is 


applied to the problem. Thus, problem (41) can be rewritten as follows: 
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s 


s 
min d. wi (cj.2) + d. wi (diy), 
= ol 


s.t. Ax + Byy < ky, 
8 


So 1 20, for t= 1, x, 58, 


=U = 
w=1 
Zo = max t t 
@)E <u < WY, where : ee Pee, 
eS epas ley eat yey vesd ae min wi(c,.v) + y w,(d;-y), 
(cE ee hid 5 wok pees EY 
dj nj’ “lj? sms’ y i=1 i=1 
055. 0G5+ C47 > Sigs Rats 72, for all a, j. = = 


s.t. Aow + Boy < ko, 
t 
Sow} =1,w, >0, for i=1,...,t. 
i=l 
(61) 
Since model (61) finds the maximum of all minimum objective values, we 
should consider the objective function of leader in the worst case, that is, 
we set c; and d; to their upper bound, in other words, (é)¥ and (di)¥, 
respectively. 


s 


: ((a.\U ((d.\¥ 
nex L~ wil (Ei)a i) 2, wil (di)e -Y); 


s.t. Ax + Byy < ky, 
8 


Sa = 1 0, for 41, 01 ~, 8; 
i=1 


=U = 

Z oy = max t t 
(@)E <v < (YU, where : rt rt 
wag seg Myre shag, min Y wi(c,.") + s w,(d;.y), 
45> bij+ Cig» Sig» hits RjQ, for all i, j. yeY = ra 


S.t. Aor + Boy < ko, 


t 
Sa La > 0, fer r= Tange 
i=1 

(62) 

To obtain an upper bound (or a maximum) for the minimum objective value, 
we should consider the values of components of the matrices Ay, By, Az, Bo, ki, 
and kz in (A1)£, (By)£, (A2)%, (Bo), (k1)¥, and (k2)Y, respectively, because 
we can find the smallest feasible region, in this case. Hence, the upper bound 
is obtained for the objective function. Doing this work, problem (19) is con- 
verted to the following problem. 
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st. (Ay)E.n + (Bi)Ly < (kU, 


s 


Yoh ee Car. 


=U i=l 
Be is ee t t 
san ae fend ate ae min S wi(c.a) + S/ wi (dy), 
for all j. yeY a Par 
st. (Ag)k.a + (Bo)k-y < (ko)¥, 


t 
Sows =1,w),>0,fori=1,...,t. 
i=1 
(63) 
We apply Theorem 2 for problem (3) and rewrite the problem as follows: 


min } | wi((@)a az) + So wi((di)e y), 
t=1 i= 
pt. (Ai )e es (Biee < (ie, 
t 
S> wi(di-y) — ((k2) 8 — (Aa)h.2)? .u = 0, 
ra i=1 ; 
= max ~ 
° @ME sus wy, ((B2)£)? u > So ui.d;, 

where v=c4,,..., ong? — 
Gi gates a, for all j. i=l 


ww; > 0, fori=1,...,8, Jj=1,..., 4, 


Z,Y,U= oO. 
(64) 
t 
Let U = {u > 0 | ((Ba)E)T.u > pee CALa Suppose that the vertex 
i=1 
points of set U are denoted by a set UY = {u',...,u"’}. Hence, we should 


solve at most r’ linear programming problems as follows: 
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=U : . 2 
LP(u*) :Z,,, = min) wi((G)a-z) + D> wilde -y), (65) 
i=1 i=1 
st. (Ay)E.2+ (Bid-y < (mY, (66) 
t t 
Yo widie-y S (kale — (A2)z.2)7 uk < ST wildy, — (67) 
i=1 i=1 
s t 
So wi = 1,50 wu; = Lagu, 2 0, for t= 16..,8 = Loe 
i=1 i=1 
Ly 20, (68) 
where k = 1,...,7’. Finding the optimal values of variables x and y for each 
a € [0,1], we can estimate values (&)".« + (di)4.y for cach a € [0,1] and 
i=1,...,s. Hence, we can obtain ZY. Having Z/ and ZY, we can easily 


compute the compromise vector Z. We will illustrate the process of obtaining 
ZE, Zand Z by a numerical example, in Section 5. 


Remark 1. As it was considered the formulation of FMOLBP in Subsection 
3.1 and the solution procedure for problem (29)—(33) in Subsection 3.2, the 
fuzzy numbers can be arbitrary according to Definition l(a) and there is 
no restriction on the fuzzy number type. The solution procedure does not 
depend on the form or type of fuzzy numbers. 


In the next section, we design an efficient algorithm to solve the problem 
of FMOLBP using the points in this section. 


4 An algorithm for resolution of the problem of 
FMOLBP as (29)-(33) 


An efficient algorithm is designed to solve problem (29)—(33). This algorithm 
produces the fuzzy compromise objective values of the problem as fuzzy com- 
promise objective values in the leader and follower and fuzzy optimal weighted 
objective value. The step-by-step algorithm is as follows: 

Step 1. Determine the weight vectors W = (wi,...,ws) and W’ = 


s t 

(wi,.-.,w,) such that So wi = 1, So wi = 1, and w; = 0,w; > 0, for 
i=1 i=1 

a. ee es ee | 

Step 2. Determine the fuzzy parameters 0 = {(v,ps(v)) | v € S(0)} of 


problem (29)—(33), where @ is one of the parameters C1;,..-, Enis C1j)+++5Enjs 
diz,...,dmi, Dag ig lags 4= Les nisi 839 => Tes ccacbs Biz; i= i Pare ee) => 


Leng Mi Ope EL DIS Dye My Cie SS Lye ig, F = Wyse se Siz, 


20 Abbasi Molai 


i=1,...,qj=1,...,m, ka, kj, i=1,...,p,j =1,...,¢, and S(@) is the 
support of v. 
Step 3. Create a-cuts of the fuzzy parameters 0 as follows: 


()a = [(8)g, (@)a] = [min{v € S(6) | pa(v) > a}, 
max{v € S(0) | pa(v) > a}}. 


Step 4. Create the following system: 


Step 5. Select some discrete points for a from [0,1] as a; = 0.1 x t for 
t=0,1,...,10 or ay = 0.25 x t for t = 0,1,2,3,4. The selection depends on 
the decision-maker or the required accuracy (it can be selected one of them 
or another partition from [0,1] for choosing the values of a ). 

Step 6. Compute the vertex points of the system of Step 4 for values a; of 
Step 5. Denote the vertex points by U2, = {ul,...,u™'}. If the system 
of Step 4 is empty for cach a € [0,1], then problem (29)—(33) has no com- 
promise solutions. The feasible solution set is empty. End. 

Step 7. Create the linear programming problem (57)—(60), for UZ, = 
fubo....,u™t} and value a;. Solve r linear programming problems and 
select the minimum objective value among them. 

Step 8. Find the lower bound of a-cut of the fuzzy compromise objective 
values of leader and follower or the fuzzy optimal weighted objective value 
using Step 7, for values a;’s of Step 5. 

Step 9. Create the following system: 


Step 10. Compute the vertex points of the system of Step 9 for values a; of 
Step 5. Denote the vertex points by (U’)?, = {(u’)h*,..., (u’)"}. Tf the 
system of Step 9 is empty for each a € [0,1], then problem (29)—(33) has no 
compromise solutions. The feasible solution set is empty. End. 

Step 11. Create the linear programming (65)—(68), for (U’)2, = {(u’)'**,... 
(u’)"%*} and value a;. Solve r’ linear programming problems and sclect the 
maximum objective value among them. 

Step 12. Find the upper bound of a-cut of the fuzzy compromise objective 
values of leader and follower or the fuzzy optimal weighted objective value 
using Step 11, for values a;’s of Step 5. 

Step 13. Create the fuzzy compromise objective values of leader and follower 
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or the optimal weighted objective value by two Steps 8 and 12. 
Step 14. End. 

Some numerical examples are presented to illustrate the algorithm and 
compare it with other existing approaches, in Sections 5 and 6. The compar- 
ison shows the efficiency of the proposed algorithm. 


5 Numerical example 
In this section, we illustrate the proposed algorithm in Section 5 to solve 
FMOLBP by a numerical example. 


Example 1. Consider the following problem: 


Z =min F(a, y) = (@.4 + dy-y, G0 + doy)’, 
rEeX 

S.t. Ay.x + Buy < ki, 
min f(x,y) = (G0 + dy, ca + dy)”, 
y 


S.t. Ao.2 + Boy < ho, 


(1,2, 3), ds = [da] = [(-1, 1, 2)], 2 = [é12] = 


= [( 
y, & = [4] = [2,3,4)], a) = [di] = [( 
)], ds = [dg] = [(1.5,2.3,3.5)], Ar = e& = 
z 
9 


] 
(ap) = 2) = (5482) = EE) = (BE). 
Ao = [éu] = [(1,1,1)], Bo = [Su] = [(2.2,3.4,4.6)], and 


Step 1. The weight vectors are wy = wo = wi = ws = 0.5. 

Step 2. The fuzzy parameters 0 = {(v, ws(v)) | v € S(@)} of problem (29)— 
(33) for this example are given above. 

Step 3. The a-cuts of the fuzzy parameters of 6 = (v1, v2,v3) can be com- 
puted the following formulas: (0) = [(#)%, (@)¥], where (6)& = v1 +(v2—-u1)a 
and (0)¥ = v3 — (v3 — v2)a. 

Step 4. The system corresponding to this step for this example is 

U ={u>o| (4.6 —1.2a).u > w{(5 — a) + w9(3.5 — 1.2a)}. 

Step 5. The selected values for a are as 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 
0.8, 0.9, and 1. 


Step 6. The vertex points of set U% are easily obtained for each a = 


O,0.1,.+0) 188 uy = 28a e000 5-120 


. Moreover U® is singleton for each 
a. 
Step 7. The linear programming problem (57)—(60), for U2, and the value 


a is created and the problem related to Z. , 1s as follows: 
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LP(u*) Zo p = min(wi(1 + a) + we(2+a))x + (wi(—1 + 2a) + we(44+ a))y, 


(-1-a) (-a+1) (5 + a) 
— ( (—a) ) oe & + -) US (‘3 + 7 
(w,(3 + a) + wh(1.5 + 0.8a))y < ((5 +a) —z)uUa 
< (w} (5 — a) + wh (3.5 — 1.2a))y, 
zy = 0. 


Step 8. The lower bound of a-cut of the fuzzy compromise objective values 
of leader using Step 7, for values a;’s of Step 5, is solved. Its results display 
in Table 1. 

Step 9. The system corresponding to this step for this example is 

U' = {u>o| (2.24+1.2a).u > wi (5 — a) + w5(3.5 — 1.2a)}. 

Step 10. The vertex points of set (U’)® are easily obtained for each a = 


0,0.1,...,l as ul, = aaa ee Moreover (U’)® is singleton for 


each a. 
Step 11. The linear programming problem (65)—(68), for (U’)?, and the 
U 


value a is created, and the problem related to Zax is as follows: 


LP(u*) : 


Zk = Min(wy.(3 — a) + w2.(5 — 2a)).2 + (wr.(2 — a) + w2.(6 — a)).Y, 
(-3+a)) | ((-1+a)) 2 (7a) 

a (4) a+ (GEE) 12-2). 

(wi (3 + a) + w)(1.5 + 0.8a))y < ((8 — 2a) — x).u), 


< (w} (5 — a) + wh (3.5 — 1.2a))y, 
Ly = 0. 


Step 12. The upper bound of a-cut of the fuzzy compromise objective values 
of leader using Step 11, for values a;’s of Step 5, is solved. Its results display 
in Table 1. 

Step 13. The fuzzy compromise objective values of leader by two Steps 8 
and 12 are created in Table 1. End. 

Figure 1 displays the fuzzy objective values of leader using the data of Table 
1. Using the data of columns 1, 2, and 3, the fuzzy objective values of follower 
can be easily computed. 


In the next section, the fuzzy compromise objectives of leader and follower, 
for an example from [3] as a benchmark, are computed and compared with 
other existing methods. It shows the efficiency of the proposed algorithm in 
this paper. 
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i : 5 a : 5 \L oe 
Table 1: Estimation of compromise vector 7 = (41, Z2) using (Zi). and (Zi). for 
i= 1,2, with different a-level values. 
a (x*,y*) optimal | (z*,y*) optimal vA Zo 
of problem re- | of problem re- 
<L =U 
lated to Z, lated to Z, 
Sake Sau Sade a 
@- ex. [ee 
(@)e r+(di) y (@)e a+(d1) y (@)¢ r+ d,) y (@)u r+(d2) y 
o a o a 

0 (0,1.087) (0,3.636) -1.087 7.272 4.348 21.816 
0.1 | (0,1.138) (0,3.362) -0.9104 6.3878 4.6658 19.8358 
0.2 | (0,1.193) (0,3.115) -0.7158 5.607 5.0106 18.067 
0.3 | (0,1.25) (0,2.891) -0.5 4.9147 5.375 16.4787 
0.4 | (0,1.311) (0,2.687) -0.2622 4.2992 5.7684 15.0472 
0.5 | (0,1.375) (0,2.5) 0 3.75 6.1875 13.75 
0.6 | (0,1.443) (0,2.329) 0.2886 3.2606 6.6378 12.5766 
0.7 | (0,1.516) (0,2.171) 0.6064 2.8223 7.1252 11.5063 
0.8 | (0,1.593) (0,2.025) 0.9558 2.43 7.6464 10.53 
0.9 | (0,1.676) (0,1.89) 1.3408 2.079 8.2124 9.639 
1.0 | (0,1.765) (0,1.765) 1.765 1.765 8.825 8.825 


T T T T T T T T T 7 T T T T T T T T T T 
Blab le Lele] lel ele Lela | ele] ENG] e@lad ela] ela] ela 
Lit i09 te eititi ti titi titetiti titi tet tt ti tt 
ie adaidadadadyaqye?yae?yedepeyaede dapat atadadadeiqyaedidia 
ent peg fe} fe 
1 1 1 1 1 1 ' 1 ' 1 1 1 1 1 1 1 i i 1 i ! 1 1 1 
itibalitigifififipifititityg ti fifififititipi ti | i 
erence eee | ele ee | eye Se sya eye ele] eee 
' ' ' ' ' ' ' 1 ' 1 ' ' ' t ' ' ' ' ' ' ' i) ' 
' ' 4d ' ' ' ' ' ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' ' ' 
BBO 3 
ti ehelalaelaielalails | ela lela g]s | eladel’ | ela] al 4 
fe fa-} a pt a fa a fe fa a at 
iyi dat adidapaepedayayayayadya dai pidayapedayeqpayadya 
itebititititiga ti: titititie titi titi tte tit ta 
ida da Papaya dedi ?didipegpeda ye dada dadapadyadedvaidada 
ideipatitititisprititititisititiftititifpipi titi ti | 
i 1 1 1 i i i 1 i i 1 i i i 1 1 i i i i ' i i i 
1 ' ' 1 ' ' ' ' 1 ! 1 ! ! ' 1 1 ! 1 ' ' ! 1 
tla pile lela tela die | cap ee la] ela] tla ea a] 
pe ee ee ee ee ee ee ee ee ee ee re ee ee ee ee ee ee ee : ZF 
2410312 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 


Figure 1: Membership function of compromise vector Z= (Z1, Z2), where Z, and Z 
are computed by ((Z)¢ ; (Z)") and ((Z2)q ; (Z2)/)), respectively. 
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6 Comparison with other existing methods 


In this section, we present an example from [3, 4] as a benchmark to show the 
efficiency of the proposed method in this paper with the existing methods. 


Example 2. Consider the following problem: 


f(z) = (x1 + 3x9 + 2ar3 + 324) 
min | fio(a) = (221 + 9x + 3x3 + 524) 
21,02 ~ ~ ~ ~ 
fis(x) = (8a1 + 9xq + 9x3 + 24) 
S.t. min (Be 2 (62 v 3x2 al 2a3 T a) 
3,L4 fro2( j= (521 - 9x9 = 9x3 + 624) 
S.t. 321 —%@2+7%3+ 324 < 48, 


201 + Ars + 2x3 _ 204 < 35, 


x1 + 2x9 — 23 +24 < 30, 


U1, 02,23, 04 2 0. 


The fuzzy numbers in [3, 4] are assumed to be triangular fuzzy numbers and 

given as follows: 

2 =4(0,2,2),8 = Oye Oe ee 

(6,8, 10),9 = (8,9, 10), 3 0 = (28 , 30, 32), 35 = (33, 35, 37), and 48 = (45, 48, 49). 
= —<U 


The problems related to Z.. and Z. are as follows: 


+ (2a) +(2+a)+ (8+ 4))23 + 3(2 +a)+ (4+ a) + 1)xa, 
s.t. 0.5((2a) + (—8 — a))x3 + 0.5((2a) + (5 + @)) 24 
< ((45 + 3a) — (4-a)a + x2)uh 
+ ((33 + 2a) — (3 — a), — (5 — a)aq) uk 


%1,%2,%3,V%4 2 0, 
where u*’s are the vertex points of the set 
U = {u = (ui, u2,u3) > o | ur + (38 — a)ug + ug > -2.5—a,(4-— a)ui 4 
(—2a)uz — uz > 5 — a} and 
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LP(u*) Z,,4, =min (1 t(3—a)+(4-—a))a1 


(4 a) + (10 — a) + (10 — a) arg 
(G3 a) + (4—a) + (10—a))as 
(Ede paella ea, 


s.t. oka 10 + 3a)r3 + (5 + 8a)ax4) 
< ((49 — a) — (2+ a)a1 + 22)uh 
((37 — 2a) — 20.21 — (3 + a) x2) us 
((—28 — 2a) + 21 + (3 — a)ax2) 
< 0.5((—5 — 2a)x3 + (10 — 2a)r4), 


U1, 02,3, 04 = 0, 


where u*’s are the vertex points of the following set: 

U = {u = (us, v2, u3) > o| uw + (Qa)ug + ug > —2.5-—a, (2 +a)ui + (-3 4 
a)uz — ug > 5— ah. 

In the above models, a € [0, 1], wi = w2 = w3 = §, and wi = wh = 0.5. The 


compromise values _ fits fi Flas fon, and fo2 are computed using the values 
of (fii)ay (fit)a (fia )a) fia)» Fta) es Fis). Gara ( fan) e » (foa)es (faa); 


and x* obtained from solving two models LP(u"). The optimal solutions 


of problems related to Ze and VAY for different values a are presented in 
Table 2. The values of leader and follower objectives for different values a 
are showed in Tables 3 and 4, respectively. 


Table 2: Optimal solutions of problems related to Zh and Zu for different values a. 


a x* optimal of prob- | x* optimal of prob- 
lem related to Ze lem related to Aan 

0 (11.25,0,0,0) (24.5,0,0,0) 

0.25 | (12.2,0,0,0) (21.66667,0,0,0) 

0.5 (13.28571,0,0,0) (19.4,0,0,0) 

0.75 | (14.53846,0,0,0) (17.54545,0,0,0) 

i (16,0,0,0) (16,0,0,0) 


A comparison among the proposed method and the methods in [3, 4] is 
presented in Table 5. 

As it is shown in Table 5, the results obtained from the proposed method 
in terms of fuzzy numbers are more accurate than methods of finding a so- 
lution with the highest membership degree. In the proposed method, the 
decision-maker can consider all of the solutions in all the membership de- 
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Table 3: Values of leader objectives for different values a. 
a Jia fiz fis 
= LD ws UT = LD rz U0 eS LD z UT 
(fu) | (fn) | (fz) | (fa) | (fis) | (As) 
0 11.25 24.5 0 73.5 22.5 98 
0.25 12.2 21.66667 6.1 59.58334 27.45 81.25001 
0.5 | 13.28571 19.4 13.28571 48.5 33.21428 67.9 
0.75 | 14.53846 | 17.54545 | 21.80769 | 39.47726 | 39.98077 | 57.02271 
1 16 16 32 32 48 48 
Table 4: Values of follower objectives for different values a. 
a far f22 
7 L ~ U = LD ry U 
(fa) (fx) (fz2) (fez) 
0 56.25 171.5 45 147 
0.25 64.05 146.25 51.85 124.5834 
0.5 | 73.07141 126.1 59.7857 106.7 
0.75 | 83.59615 | 109.6591 | 69.05769 | 92.11361 
1 96 96 80 80 
? 
! | 
+ {\ + 
, [| . ‘i lL \ 
j ‘| kl fp \ i r 
| \ | \ j \ 
4 f \ ' ‘ 
| \ | \ | 
pL ‘ , i—\ 


Figure 2: Fuzzy compromise objective values of leaders by the proposed approach 


Table 5: Comparison among the proposed method and two other methods. 


The method of The method of The proposed method 
Baky Eid, and El Pramanik and Dey 
Sayed, in [3] in [4] 
fu = 29 fy =1 fu = 37 My = 0.902 fir in Figure 2 

fig = 70.9483 | py2 = 0.917 fie = 90 baz = 0.815 fig in Figure 2 

fig = 88.2743 | fiz = 0.821 | fig = 108.25 | pig = 0.692 fig in Figure 2 

for = 80.912 | proy = 0.527 | fo, = 78.25 | po, = 0.496 foi in Figure 3 

foo = 111.27 | p22 = 0.81 | foo = 105.5 | poo = 0.795 fog in Figure 3 
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Figure 3: Fuzzy compromise objective values of followers by the proposed approach 


grees. Moreover, the method of Pramanik and Dey in [4] did not follow the 
basic concepts of bilevel programming. Also, the obtained compromise solu- 
tions improve the compromise solutions in [3] for degrees 1, 0.917, 0.821, and 
0.81, considerably. Moreover, the obtained compromise solution covers the 
compromise solution in [3] for degree 0.578. 


7 Conclusion 


This paper showed how Zadeh’s extension principle can be efficiently applied 
to solve the FMOLBP problem. Using the structure of bilevel programming 
and concept of extension principle, two crisp multiobjective linear three-level 
programming problems were designed to compute the lower and upper bound 
of the a-level of the compromise objective value. Due to the strongly NP- 
hardness and high computational complexity of the problem, there is no 
efficient algorithm to solve it. The developed methods cannot completely 
solve the problem. To overcome the deficiencies, a pair of crisp multiobjec- 
tive linear three-level programming problems was proposed to construct the 
compromise fuzzy objective functions of leader and follower. The main ad- 
vantages of this method are access to fuzzy compromise objectives of leader 
and follower with less computational complexity with respect to approaches 
based on the rejection of solution and resolve model successively until find- 
ing a satisfactory decision or based on FGP. The proposed method is based 
on the weighting method, dual theory, and vertex enumeration approaches 
that preserve the linearity of the model. A comparison among the proposed 
method and the methods of Pramanik and Dey [4] and Baky Eid, and El 
Sayed, [3] showed the efficiency of our method. 
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